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Abstract. It has been recently suggested that magnet- 
ically affected neutrino oscillations inside a cooling pro- 
toneutron star, created in a supernova explosion, could 
explain the large proper motion of pulsars. We investigate 
whether this hypothesis is in agreement with the observed 
properties of pulsars and find that present data disfavor 
the suggested mechanism. The relevance of our results for 
other models proposed to understand the origin of pulsar 
velocities is also discussed. 
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1. Introduction 

One of the challenging problems in pulsar astrophysics is 
to find a consistent and observationally verified explana- 
tion for the high peculiar velocities of pulsars. These veloc- 
ities can be as high as 1000 km/s and have a mean value of 
450 km/s, much greater than the random velocities of or- 
dinary stars (Harrison et al. 1993| ; Lyne & Lorimer 1994 ). 
Several mechanisms have been put forward in the past to 
explain the origin of the large proper motions. Since it 
is believed that pulsars are born during the first stages 
of some Type II or core-collapsing supernovae, an asym- 
metric core collapse or explosion could give the pulsars 



the observed velocities (Shklovskii 197C|; Wo osley 1987; 
Woosley & Weaver |l99|; Janka & Miluer |1994 Burrows & 
Hayes 1996). The evolution of close binary systems could 



also be responsible for the large pulsar velocities (Gott et 
al. 1970). Alternatively, the emission of electromagnetic 



radiation during the first months after the supernova ex- 
plosion, stemming from an off-centered rotating magnetic 
dipole in the newborn pulsar, could g ive th e pulsar a sub- 
stantial kick (Harrison & Tademaru 1975 ). Another ap- 
proach is based on the assumption that most of the en- 
ergy and momentum released during a Type II supernova 
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explosion (~ 10^^ erg) are carried off by the relativistic 
neutrinos, as was observationally confirmed by the detec- 
tion o f a ne utrino burst a ssocia ted with SN1987A (Hirata 
et al. 1987 ; Bionta et al. 1987 ). Therefore, an asymmet- 



ric neutri no em ission, caused for e xamp le by convection 
(Woo sley |l987t Woosley & Weaver |l992[ Janka & Midler 
1994 ) or s trong ma gnetic fields (Chugai 1984; D orofee v 
et al. |1985| ; Vilcnkin |l995| ; Horowitz & Piekarewicz |l997|) , 
may play an important role when trying to understand 
the origin of pulsar velocities. Not all these mechanisms, 
however, seem to be able to produce the required high ve- 
locities and many of them suffer from the fact that they 
need magnetic fields of the order of > 10^^ G, pulsar pe- 
riods of ^ 1 ms or that they are not sustained by other 
observed pulsar properties (Duncan & Thompson 1992). 

In a recent paper Kusenko & Segre ( |l996a[ ), hereafter 
KS, have proposed a new mechanism of asymmetric neu- 
trino emission based on resonant neutrino oscillations. The 
three types of neutrinos, i^fj, and Vt, are abundantly 
produced in the core of a collapsing star which later on 
may become a pulsar. The matter density is so high in 
the core that the neutrinos do not escape but get trapped. 
They undergo a diffusion process until they reach higher 
radii, where the density has decreased and they can freely 
stream away. The emission surface, the so-called neutrino 
sphere, is not the same for the three types of neutrinos. 
Since electron neutrinos can interact via both charged 
and neutral currents they interact more strongly in the 
protoneutron star than muon and tau neutrinos. Hence, 
the electron neutrino sphere is at a larger radius than 
the muon and tau neutrino spheres. The authors in KS 
showed that under these conditions neutrinos i^r can res- 
onantly turn into v^, by means of the adiabatic MSW 



effect (Smirnov 1996), in the regioiL. between the tauonic 
and the electronic neutrino spheresEJ. The emerging elec- 
tron neutrino, however, will be absorbed by the medium 
and therefore the resonant surface becomes the effective 
surface of emission of the Vr- Neutrinos propagating in 



For definiteness, only the two flavors Ue and Ut have been 
discussed with Am^ = 771^(1/^) — rn^ive) ~ rn^ivr) and small 
mixing. 
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media with a longitudinal magnetic field B have differ- 
ent electromagnetic properties than in the vacuum case. 
They acquire an effective electromagnetic vertex which is 
induced by weak interactions with the charged particles 
in the background and generates a contribution oc B ■ k 
to the effective self-energy of the neutrin o, k b eing the 
neutrino mome ntum (Esposit o fc C apone 1996 ; D'Olivo 
Elmfors et al. 



& Nieves 1996 



1996). The induced vertex 



modifies the flavor transformation whilst preserving chi- 
rality and, as a result, the location at which the resonance 
occurs is affected, leading to the spherical symmetry of 
the effective emission surface being replaced by a dipolar 
asymmetry. The condition for resonant oscillations to take 
place is accordingly given by 



2k 



■cos 261 = V2GFNe{r) + 

eGp f 3N^{r )Y^^ k- B 



(1) 



where 9 is the neutrino vacuum mixing angle, Gp the 
Fermi constant, iVc(r) the charge density of the degenerate 
electron gas in which the neutrino propagates and r the 
radial coordinate. Neutrinos emitted from the two mag- 
netic poles of the resonant surface then have slightly dif- 
ferent temperatures because the two poles are at slightly 
different radii. The outcome is an asymmetric emission of 
momentum carried by neutrinos which gives the neutron 
star a kick in the direction of the magnetic field and thus 
leads to a recoil velocity in agreement with observational 
data. Quantitatively the kick is described by the asymme- 
try in the third component of momentum and estimated 

by 



M /3MeV 
k V T 



B 



3 X IQi^G 



(2) 



Since the total momentum carried away by neutrinos emit- 
ted by the protoneutron star is 100 times the momen- 
tum of the proper motion of the pulsar, an asymmetry 
of 1% would give a kick in agreement with observation. 
Assuming an average energy for the tau neutrinos leaving 
the protoneutron star of « 10 MeV, which corresponds to 
T w 3 MeV, the authors in KS obtain the desired asym- 
metry of 1% for values of i? ~ 3 x 10^^ G. As an ad- 
vantage over other neutrino emission mechanisms the one 
discussed here works for smaller magnitudes of the mag- 
netic field and does not demand for any constraints on the 
pulsar period. If the resonant neutrino conversion turned 
out to be the origin of pulsar velocities, one could use pul- 
sar observations to obtain information on neutrino masses. 
The implications for particle physics models would be of 
great importance. 

A possible weak point of the model, however, lies in 
the fact that the mass of the tau neutrino is required to 
be rrii,^ ^ 100 eV in order to have a resonant conversion 
in the protoneutron star between the electron and the tau 



neutrino spheres. Although such a mass is not excluded 
by laboratory experiments, it is difficult to accommodate 
it in the standard Big Bang cosmological model. Indeed, 
the present age of the universe can be used to set an up- 
per bound on the mass of any light, stable neutrino species 



(Gerstein & Zel'dovich |1966| ; Cowsik & McClelland [197^ : 
< 92 eV, where h is the Hubble constant in units 
of 100 km s^^/Mpc. Nevertheless, this problematic point 
might be overruled by unstable tau neutrinos or by a bet- 
ter understanding of the evolution of the universe. An- 
other comment can be made in the context of the super- 
nova physics. From treatments of neutrino transport it 
was found that the average energy for tau neutrinos leav- 
ing the supernova can be as high as 27 MeV (Janka 1992). 
Such a value, however, would demand for a higher mag- 
netic field in Eq. (^) in order to still explain the observed 
pulsar velocities thus weakening one of the advantages of 
the proposed model. On the other hand, it is interesting to 
note that neutrino oscillations between Vg and a tau neu- 
trino with a mass of 10— 100 MeV could, for a large range 
of mixing angles, help to explode supernovae while leaving 



r-process nucleosynthesis unscathed (Raffclt 1996). A bet- 
ter understanding of supernova explosions and data from 
future galactic supernovae could be used to test whether 
the asymmetric resonant neutrino emission model is actu- 
ally realistic. 

Since none of the above mentioned points represent a 
serious danger to the discussed mechanism, but its impli- 
cations for astrophysics and physics beyond the standard 
model of the fundamental interactions would be signifi- 
cant, it is worthwhile to examine further related observa- 
tional consequences with the aim of arriving at indepen- 
dent constraints. In addition, the derived constraints can 
be used to test the viability of some of the other proposed 
models for the origin of pulsar velocities. In this paper 
we derive, using the mechanism introduced in KS, a new 
expression relating the transverse velocity, the magnetic 
field and geometric parameters of fast spinning pulsars. 
We analyze observational data on pulsars to see whether 
this relation is fulfilled. Our work, however, seems to indi- 
cate that pulsars do not satisfy this relation, and therefore 
present data do not support the neutrino oscillation mech- 
anism. 

2. Observational tests 

According to the authors of KS the asymmetric neutrino 
emission from a magnetically distorted neutrino sphere 
gives a net momentum to a pulsar along the magnetic 
axis. It is believed that, in general, the spin and magnetic 
axes of pulsars are not aligned, and thus the magnetic field 
rotates around the spin axis with the period P of the pul- 
sar. From this it can be seen that a net velocity along the 
magnetic axis is only obtained if the pulsar period is much 
larger than the characteristic time scale on which the 
neutrino flux from the cooling protoneutron star remains 
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constant. For such a case, the length of the velocity vector 
V and magnetic field vector B in the core of the pulsar 
should be proportional: 



the two groups and, for calculating the angles, they find 
different empirical laws relating the pulse width and the 



V (X B . 



(3) 



In a different paper Kusenko & Segre (1996b) apply the 
model to pulsars which are slow spinning at birth, P > 1 s 
and P > 1.3 s, and argue that this B—v correlation should 
become increasingly more significant as the period of ro- 
tation approaches a few seconds. They claim that observa- 
tional data show such an increase in the correlation for the 
examined periods, which would support the investigated 
model. 

In the present paper we consider fast spinning pulsars, 
P < 1 s, instead. The characteristic time t,y is of the or- 
der 1 — 10 s (Suzuki 1994). Therefore, a reliable test of the 



mentioned correlation between B and v for slow spinning 
pulsars should only be possible for periods P ^ 10 s which 
are not observed. We believe that pulsars with P < 1 s 
form a better set to test the model. A period P < 1 s is 
sufficiently smaller than the time interval in which neutri- 
nos are emitted and therefore it can be safely assumed that 
the magnetic field and thus the net momentum vector ro- 
tate around the spin axis several times during the neutrino 
emission. The momentum component on the plane orthog- 
onal to the spin axis averages out and only the component 
parallel to the spin axis survives. Thus, for fast spinning 
pulsars, the proper motion should be in the direction along 
the pulsar spin axis. As an important consequence of this 
alignment of the velocity and spin axis for fast spinning 
pulsars one can easily prove that the relation 



vt oc B cos a sin(a -t- (3) 



(4) 



must apply where vt is the projected pulsar velocity in the 
sky, a the angle between the magnetic and spin axes, and 
P is the angle between the magnetic axis and the line of 
sight (impact parameter). We will see that observational 
data allow a test of this relation. 

Pulsar proper motions are determined by radio inter- 
ferometry (Harrison et al. 1993|) a nd by interstellar scintil- 
lation observations (Cordes 1986| ). In order to determine 
the emission geometry parameters a and (3 one has to 
rely on some emission model for pulsars. The so-called 
magnetic-pole model assumes that the electromagnetic ra- 
diation originates in the vicinity of a magnetic pole. The 
radiation is emitted in a narrow cone-shaped beam cen- 
tered around the magnetic axis. The nonvanishing an- 
gle a between the spin axis and the magnetic axis pro- 
duces the characteristic pulses every time t he be am sweeps 
the line of sight (Lyne & Graham-S mith 1990 ). Tw o dif- 
ferent groups (Lyne & Manchester |l988| ; Rankin |l990 , 
1993a| , [I993b|) have determined a and /3 for more than one 



hundred pulsars. The method used by both groups is based 
on the apparent increase of pulse width as 1/sina. The 
classification scheme of pulsars, however, differs between 



pulsar period: while Lyne & Manchester (1988) use width 
cx p-i/3, Rankin ( [1990| , [19934 |1993b| ) employs width 
cx P~^/^. The results of both groups show, in general, 
a remarkable agreement for angles a < 40° but t hey c an 
widely differ for larger ang les (Miller & Hamilton [l99^ . 

Our aim has been to find out whether the relation 
given by Eq. (^) is corroborated by observational data. 
Bearing in mind the lack of complete agreement between 
the different measurements of the emission angles, we have 
studied three different sets of data: (a) angles f rom Ly ne & 
Manchester ( 1988 ), (b) angles from Rankin ( 1993b| ) and 
(c) the average of the angles from sets (a) and (b) for 
those pulsars which are common in both of these sets and 
which have similar a values, namely Aa < 12°. Further- 
more, it has to be taken into account that Eq. (^), true at 
the pulsar birth, may be spoiled by the temporal evolu- 
tion of pulsar parameters (angles and magnetic field). It 
is believed that the magnetic field of a pulsar decays with 
a characteristic time ~ 10^ years (Harrison et al. 1993). 
There is also observational evidence supporting that the 
magnetic and spin axes become aligned when the pulsar 
ages (Lyne & Manchester 1988 ). To avoid these problems 
associated with time evolution we have selected from the 
data the pulsars younger than 10^ years and, to be even 



safer, younger than 3x 10^ years (Lorimer et al. 1995) thus 
creating two subsets for each of the sets (a), (b) and (c). 
We disregard older pulsars since their present properties 
may be rather different from those at birth. We do not 
consider pulsars in binary systems either since for them 
Eq. (^ could be affected by the evolution of close binary 
systems. The selected pulsars for the three cases (a), (b) 
and (c) are listed in Tables |l| and |. 



Table 1. Pulsars of case (a) as described in the text. The 
pulsars with an asterisk have ages of less than 3 Myr, all the 
others of less than 10^ yr. The pulsars which are also members 
of set (c) are primed. 



PSR B0136-I-57*' 
PSR B0458-f46" 
PSR B0736-40' 
PSR B1642-03' 
PSR B1933-I-16*' 



PSR B0329+54' 
PSR B0559-05 
PSR B0823+26 
PSR B1706-16*' 
PSR B1946+35*' 



PSR B0450-18*' 
PSR B0656-(-14* 
PSR B1508-f55* 
PSR B1929-(-10 
PSR B2020-H28*' 



We plot log [vt/ cos Of sin(a -I- /?)] versus logP for 
young, fast spinning pulsars (P < 1 s) in Figs. |], § and |[ 
for the cases (a), (b) and (c), respectively. The transverse 
velocity vt is calculated from the measured angular ve- 
locity using the pulsar distance D derived from its 
dispersion measure. The pulsar velocities have been cor- 
rected for galactic rotation by removing a flat rotation 
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Table 2. Pulsars of case (b) as described in the text. The 
pulsars with an asterisk have ages of less than 3 Myr, all the 
others of less than 10^ yr. The pulsars which are also members 
of set (c) are primed. 



PSR B0136+57*' 
PSR B0450-18*' 
PSR B0736-40' 
PSR B0919+06* 
PSR B1642-03' 
PSR B1822-09* 
PSR B1946+35*' 
PSR B2217+47 



PSR B0329+54' 
PSR B0450+55* 
PSR B0740-28* 
PSR B1449-64* 
PSR B1706-16*' 
PSR B1842+14 
PSR B2020+28*' 
PSR B2224+65* 



PSR B0355+54* 
PSR B0458+46*' 
PSR B0823+26 
PSR B1508+55* 
PSR B1818-04* 
PSR B1933+16*' 
PSR B2021+51* 



curve of rotational velocity 225 km/s, with the sun at a 
distance of 8.5 kpc from the galactic center. They have 
also been corrected for the peculiar solar motion of the 
sun which was assumed to to have a velocity of 15.6 km/s 
in the direction I = 48?8 and b = 26°3 (Murray |l98|). 
The parameters D, B, P and the pulsar ages as well as 
the parameters needed to compute fix have been taken 
from Taylor et al. (199E, 1995), an updated catalog con- 
taining 706 pulsars. 
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Fig. 2. Same plot as Fig. |l| but now the angles a and f3 taken 
correspond to case (b). 
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Fig. 3. Same plot as Fig. ^ but now the angles a and /3 taken 
correspond to case (c). 



Fig. 1. This figure shows log [tiT/(cos a sin(Q + /3))], where vt 
is in km/s, versus the surface magnetic field logB, where B is 
in Gauss, for pulsars with period P < 1 s and age less than 
10^ yr (all circles). Those pulsars among them which have an 
age of less than 3 x 10® yr are indicated by filled circles. The 
angles a and /3 taken correspond to case (a). The logarithm is 
to base 10. 



The scatter of points in the plots we depict suggests 
that there is hardly any correlation. To evaluate rigor- 



ously their significance, we use the Spearman rank-order 
correlation coefficient rg and its probability 'P{rs)- The 
coefficient rs ranges from 1 (perfect correlation) to — 1 
(perfect anticorrelation) . Further information is given by 
V{rs), the probability that rs for two uncorrelated data 
sets would be larger than the Spearman coefficient found. 
Hence, to have a significant correlation between two sets, 
one must obtain rs ^ I and V{rs) <C 1. The calculated 
rg and Virs) corresponding to our three cases are shown 
in Table |[ We calculate the two numbers using pulsars 
younger than 3 x 10® yr and pulsars younger than 10^ yr. 
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We are unable to find any correlation of the sort given 
by Eq. (except possibly for the case (b) with pulsars 
younger than 3 x 10^ yr). 



Table 3. Spearman rank-order correlation coefficient rs and 
associated probability V{rs) for the three cases discussed in the 
text, considering pulsars younger than 10^ yr and 3 x 10® yr. 
The number of pulsars in each set is also given. 



Case 


Age[10^ yr]< 


Number 


rs 


V{rs) 


(a) 


1 


15 


6.1 X 10"^ 


0.83 


(a) 


0.3 


9 


0.25 


0.52 


(b) 


1 


23 


0.31 


0.15 


(b) 


0.3 


17 


0.68 


2.9 X 10"^ 


(c) 


1 


10 


-0.10 


0.78 


(c) 


0.3 


7 


7.1 X 10"^ 


0.88 


It 


is interesting to point out the existence of fast spin- 



ning pulsars like PSR B0833-45 which have large peculiar 
velocities whilst the magnetic and spin axes are nearly 
perpendicular. Clearly, the mechanism presented in KS 
cannot account for the proper motion of these pulsars. We 
have also examined the case of young, fast spinning pulsars 
with small angles a < 15°. In the framework of the stud- 
ied mechanism those pulsars should have small transverse 
velocities vt (since (3 is small for most pulsars). For the 
mean transverse velocity of the pulsars fulfilling the selec- 
tion criteria we obtain 344 km/s. This has to be compared 
with the mean average transverse velocity for pulsars with- 
out any restrictions on a or the pulsar period. Such mean 
values were computed to be 345 ± 70 km/s for 29 pulsars 
which are younger than 3 Myr, and 300 ± 30 km/s for 9 9 
pulsars without the age constraint (Lyne & Lorimer |1994 ) 



The fact that there is no significant difference between the 
mean transverse velocity we obtain and the latter ones 
above does not support the model examined here. 

There is yet another important consequence of the 
alignment of the spin axis and the velocity vector for fast 
spinning pulsars, which has been explored some time ago 



(Tademaru 1977; Anderson & Lyne 1983) in the context 
of the asymmetric emission of electromagnetic radiation 
mentioned in the introduction of this paper. If these two 
vectors are aligned, the difference between the projected 
pulsar velocity angle in the sky ij}^ and the position an- 
gle of the linear polarization vector in the center of a 
pulse (corrected for Faraday rotation) Vp should be ei- 
ther 0° or 90° depending on the emission mechanism and 
the inclination of the spi n and magnetic axes to the line of 



sight. While Tademaru ( |1977D found that the histogram of 
\ipv — V'pl showed two peaks at 0° and 9 0°, a subsequent 



study performed by Anderson & Lyne ( 1983 ) concluded 
that such peaks were not present. The latter result was 
argued to be in disagreement with the mechanism intro- 



duced by Harrison & Tademaru (1975). It would not sup- 
port the resonant neutrino emission model discussed here 
either. Given the significant increase in the amount and 
quality of available pulsar data over the last years it would 
be of interest to carry out the — '4'p\ test again. This 
could not only help to resolve the contradicting conclu- 
sions of the previous studies but also yield an additional 
test of the neutrino emission model if only fast spinning 
pulsars are selected from the updated and enlarged sets 
of observational data. Moreover, both tests have different 
sources of errors. While the — i/jpl test does not depend 
on the angles a and (3 the derivation of which is model 
dependent, it relies on the measurement of the rotation 
measure RM of a pulsar. The contribution from Faraday 
rotation, RMx A^, has to be subtracted from the measured 
polarization angle to obtain the intrinsic polarization an- 
gle 1pp. Our test relies on the measurement of a and (3 but 
does not need the value of RM for each pulsar which can 
be a significant source of errors. In this sense, both tests 
constitute complementary ways of checking whether the 
model given in KS and related models are realistic or not. 

3. Discussion and conclusions 

We have looked into possible ways of testing the new 
mechanism suggested by KS to explain the large proper 
motion of pulsars. 

It was argued that for fast spinning pulsars this mech- 
anism forces the vector velocity of a pulsar to be aligned 
with its spin axis. Making use of this alignment we have 
found Eq. (0), which is a relation between the magnetic 
field of a pulsar, its proper motion, the angle a between 
the magnetic and spin axes and the impact parameter (3. 
This relation should be true as long as pulsar velocities 
are produced by asymmetric emission from a magneti- 
cally distorted resonant neutrino-sphere. We have stud- 
ied whether pulsar data actually show this correlation. To 
carry out this study we have had to face the fact that there 
is some controversy about the true values of a and j3 in 
the sense that the angles from two different groups show 
discrepancies for a > 40°. Therefore, we have in a first 
step considered separately the data from each group. In 
addition, we have also selected the pulsars studied by both 
groups for which there is a reasonable agreement between 
the angles measured by either group. We believe that the 
results we derive in the latter case are the most reliable. 
In this case, observational data do not show any signifi- 
cant correlation and therefore, the mechanism introduced 
in KS is unlikely to explain the proper motion of pulsars. 
We have indicated that the temporal evolution of pulsars 
may spoil any initial correlation. To diminish this effect 
we have only considered young pulsars. Another possible 
caveat is the fact that Eq. (Q) involves the magnetic field 
in the core of the pulsar whilst what is measured is the sur- 
face magnetic field. It is, however, plausible and accepted 
to assume that both fields are proportional (Manchester 
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& Taylor |1977| ; Ruderman [l99lD . Looking at Table | one 
can see that the number of pulsars with available data 
which satisfy all the requirements regarding period, age 
and availability of emission angles as well as proper mo- 
tions is small in all the cases considered. Although our re- 
sults cast doubt on the resonant neutrino emission as the 
mechanism responsible for the proper motion of pulsars, 
clearly one needs more data, and clarify the uncertainties 
related to the true emission angles, to rule out definitely 
or support this mechanism. 

As a final remark we would like to point out that our 
conclusions also apply to other mechanisms proposed to 
explain the large velocities of pulsars. This is the case for 
fast spinning pulsars in all models where the momentum 
kick is coUinear to the magnetic axis and the resulting 
velocity proportional to B. Under these conditions veloci- 
ties along the spin axis are predicted. Hence, the work by 
Chugai (|1984[) and Dorofeev et al. (|1985|) is affected. They 



study how the polarization of electrons and positrons by 
a high magnetic field in a newborn neutron star will gen- 
erate an anisotropic neutrino emission. Also the mecha- 



nism suggested by Vilenkin (1995) is disfavored for a uni- 
form magnetic field, as is the mechanism in Horowitz & 



Piekarewicz (1997). Pulsar velocities from resonant spin- 



flavor precession of neutrinos are studied by Akhmedov et 
al. ( p^97| ). A gain, this mechanism is put under pressure 
by our work as well as the model suggested by Kusenko & 
Segre ( 1997 ) where a resonant conversion to sterile neu- 
trinos induced by neutral currents is considered. 
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